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A scenario is investigated in which the leading- twist pion distribution amplitude <Ptt(x) is approx- 
imated by the pion decay constant f n for all essential values of the light-cone fraction x. A model 
for the light-front wave function ty(x,k±) is proposed that produces such a distribution amplitude 
and has a rapidly decreasing (exponential for definiteness) dependence on the light-front energy 
combination k"x/x{l — x). It is shown that this model easily reproduces the fit of recent large-Q 2 
BaBar data on the photon-pion transition form factor. Some aspects of scenario with flat pion 
distribution amplitude are discussed. 
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I. INTRODUCTION 

The pion distribution amplitude (DA) (f w (x) [H,0] is & n important function accumulating information about momen- 
tum sharing between the quarks of the pion when the latter is in its valence qq configuration. It is an inherent element of 
perturbative QCD calculations of hard exclusive reactions involving the pion. From the solution [1, H|, @ of the evolu- 
tion equation for the pion DA it follows that independently of its shape at low normalization point jiQ < 1 GeV, at large 
values of the probing momentum the pion DA acquires universal asymptotic form: ip w (x,[i — > oo) — > 6/ w a;(l — x) [5(. 
However, in practical calculations, it is very important to know what is the shape of the pion DA at moderate and low 
scales fi. The standard measure of the width of the pion DA is the value (£ 2 ) of its second moment with respect to the 
relative momentum fraction variable £ = x — (1 — x). QCD sum rule calculations [6( give large (£ 2 ) > 0.4 values for 
this moment (compared to (£ 2 ) = 0.2 for the asymptotic DA) which indicates that the pion DA is a wide function for 
/i 2 < 1 GeV 2 . Recent lattice calculations @, B give (£ 2 ) > 0.3 for \j? values in this region. A direct calculation of the 
pion DA in the Nambu-Jona-Lasinio model [§] (see also Ref. [l(|) produces the result that ipn{%) = f-n for all values 
of the momentum fraction x, i.e., that pion DA is constant. The same result was obtained in the "spectral" quark 
model [ll|. The value of (£ 2 ) for this "flat" DA is 1/3, which is compatible with the results of the lattice estimates, 
though smaller than the result of QCD sum rules. It should be noted that the usual procedure of reconstructing pion 
DA from its moments in the CZ approach (followed by essentially all other groups) is to build it as a sum of the lowest 
(two or three) Gegenbauer polynomials corresponding to multiplicatively renormalizable components of the pion DA 
evolution decomposition. Since these components have x(l — x) as an overall factor, such a procedure excludes flat 
DAs from possible models. However, this restriction on the pion DA model building is just an assumption. In the 
present paper, our goal is to analyze the photon-pion transition form factor in a scenario with flat pion DA. The curve 
for the form factor which we obtain is in a very good agreement with recent BaBar data [12| | , which is basically the 
main motivation for our investigation. 

The paper is organized in the following way. In Section II, we give an overview of the basic facts about the pion 
distribution amplitude: its definition, evolution and results concerning its shape. Section III is devoted to the study 
of the photon-pion transition form factor. We briefly describe pQCD results for this form factor, and then calculate it 
within the light-front formalism using a model wave function ^{x, k±) that reproduces flat pion DA after integration 
over quark transverse momentum k± and rapidly (exponentially) decreases for large values of the standard light-front 
energy combination k\jx{\ — x). The k± width parameter a of this wave function can be easily adjusted to produce 
a curve practically coinciding with the data fitting curve given in Ref. [HJ . This value of a corresponds to the value 
(kj_) = (420 MeV) 2 for the average transverse momentum squared, which has the magnitude that one would expect 
for the valence qq Fock component of the pion light-front wave function. We analyze the structure of the one-loop 
corrections for a flat DA in pQCD, and find out that the optimal value /i 2 = aQ 2 of the normalization scale for a flat 
DA is very small. We argue that this is an evidence that the flat pion DA should not be evolved in our calculation 
of the photon-pion transition form factor. Finally, we discuss some aspects of the flat pion DA scenario and then we 
summarize the paper. 
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II. PION DISTRIBUTION AMPLITUDE: BASICS 



A. Definition and Evolution 

The pion distribution amplitude <Ptt( x ) mav be introduced [H as a function whose x n moments 

fn = [ X n <p n (x)dx (1) 



are given by reduced matrix elements of twist- 2 local operators 

i n+1 <0| 1(0)75 {j v D Ul . . . D„Ju(0)\tt + ,P) = {P V P V1 ...P Vn }f n , (2) 

or as the fc^-integral 



ip^(x,n) = [ ^(x,k±)d 2 k± (3) 

_ " 'fei<M 2 



(2tt) 3 



of the light-front wave function ^(x, k±). The zeroth moment of y w (x) corresponds to matrix clement of the axial 
current, and is given by the pion decay constant 

l 

<p^{x) dx = f n , (4) 

which is known experimentally. In the conventions that we use, f n ~ 130 MeV. Eq.(j4| gives an important constraint 
on the pion distribution amplitude (DA), fixing the integral under the ip v (x) curve, but it puts no restrictions on its 
shape. In fact, the pion DA depends on the renormalization scale \x that is used to define matrix elements of twist-2 
local operators: <pn(x) — > ip v (x, fi). The evolution equation for the pion DA may be written either in matrix form [l[ 

d " 

M"T" /nO) = Z nkfk(l*) (5) 

" k=0 

(see also [lH) or in kernel form Q 

The solution of the evolution equation was obtained [H, 0| , m the form of expansion over Gegenbauer polynomials 

■ f2n//3o 1 



fj,^(p^(x,fi) = I V(x,y)(p„(y,fj,)dy . (6) 



o„(x,n) = 6/,r x(l - x) { 1 + ^ a2nCl^(2x - 1) 

n=l 



(7) 



where 72n > is the anomalous dimension of the composite operator with 2n derivatives, and f3o is the lowest 
coefficient of the QCD /3-function. As a result, when the normalization scale /i tends to infinity, the pion DA acquires 
a simple form [5| 

<p„(x,fi oo) = 6/^(1 - x) , (8) 

known as the "asymptotic DA" . 



B. Shape 

The question, however, is what is the shape of the pion DA at low normalization scales fi < 1 GeV. Some qualitative 
(and maybe overly simplistic by today's standards) argumentation about a possible shape of the pion DA was given in 
our 1980 papers [HI, [HI]. Namely, in case of a system of two equal-mass non-interacting particles, ip(x) = S(x — 1/2). 
When the interaction is switched on, the DA broadens. The width T of tp(x) may be estimated as T ~ Ei n t/m q - 
Hence, for heavy mesons (e.g., for T particles), ip(x) is rather narrow since m q ^> M ~ Aqcd ~ 300 MeV, where M 



3 



is a parameter characterizing interaction strength. On the other hand, taking m Ut d ^ lOMeV for the quarks in the 
pion, we conclude that l -p- K {x) is very broad. Assuming a simple exponential model 



^{x, k±) ~ exp 



k\ + m 2 
M 2 x(l -x) 



for the light-front wave function gives 



(pn{x, fj, ~ M) = /„■ exp 



M 2 x{l - x) 



(9) 



(10) 



In this case, the pion DA <p^{x) is close to f n everywhere outside the regions < x < m 2 /M 2 ~ 10 3 and 
< 1 — x < m 2 /M 2 ~ 10~ 3 . In these regions, ip n (x) vanishes rapidly. 




FIG. 1: One-gluon-exchange diagram for pion electromagnetic form factor in perturbative QCD. 
Initially, the pion DA appeared in the perturbative QCD expression [l[ 

for the asymptotics of the pion form factor calculated through a one-gluon-exchange diagram (see Fig[T]). Here, 
xyQ 2 is the virtuality of the exchanged gluon. If one takes the flat pion DA <p w (x) = f w , both integrals in x and y 
logarithmically diverge, which means that pQCD factorization fails in this case. Evidently, the finite size R ~ 1/M of 
the pion should provide a cut-off for the x, y integral, which suggests that xyQ 2 should be substituted by something 
like xyQ 2 + 0(M 2 ), with the additional term having a meaning of the average of the squared transverse momentum 
of the quarks. One may also treat 0{M 2 ) as an effective gluon mass squared. Then the x, y integrals are convergent, 
but the integral is dominated by the region where the nominal gluon virtuality is 0(M 2 ). This means, first, that it is 
small, and what is even more important, that it is not growing with Q 2 . For these reasons, the gluon-exchange line 
cannot be treated as a part of a perturbative short-distance subprocess, in which virtualities of all the lines should be 
parametrically 0(Q 2 ). Hence, it should be absorbed into the nonperturbative part of the diagram, i.e. into the soft 
pion wave function. The pion form factor must be then calculated in some nonperturbative way. Such a calculation 
was accomplished within the QCD sum rule approach [l6l . [l7| , with the results close to experimental data. 

The same QCD sum rule approach was used [6j to calculate £ 2 and £ 4 moments of the pion distribution amplitude 
(/>„(£) (which is the original DA ^(x) written as a function of the relative variable £ = x — (1 — x) and divided by 
/„■). The value of (£ 2 ) is a quantitative measure of the width of the distribution amplitude. In particular, (£ 2 ) is 
zero for the infinitely narrow DA ip v (x) = f v S(x — 1/2), it equals to 1/5 for the asymptotic DA ([5]) and to 1/3 for 
the flat (Ptt(x) = f„ DA. The calculation of Chernyak and Zhitnitsky (CZ) [f| gave the result larger than 1/3, namely 
(£ 2 ) = 0.40 for the "bare" value that was attributed to the normalization point /z 2 = 1.5 GeV 2 and then renormalized 
to the reference scale /j, 2 = 0.25 GeV 2 , which resulted in (£ 2 ) = 0.46. Without touching a subtle point whether a 
perturbative evolution to such a low scale is justified, we can say that the CZ results clearly indicate that the pion 
DA is a wide function, and a generalized flat DA of ^(C) — a + 3(1 — a)£ 2 type could have been used as a model 
fitting the values of (£ 2 ) and (£ 4 ) obtained from the CZ calculation. However, the fitting model 

= " e 2 (i-a (12) 

was constructed from the sum of two first terms of the Gegenbauer expansion ([7]), which has x(l — x) (or 1 — £ 2 ) as 
an overall factor, thus excluding all models with DA's that do not linearly vanish at the end-points. 
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An implicit assumption of the CZ calculation is that it is sufficient to take into account only the two lowest 
condensates (GG) and a s (qq} 2 in the operator product expansion (OPE) for the relevant two-point correlator. An 
alternative attitude [l8[ is that the quark condensate (q(O)q(O)) is just the first term in Taylor expansion of the 
nonlocal condensate (q(0)q(z)) = (qq) f{z 2 ) that explicitly appears at the initial steps of OPE calculations. Modeling 
f(z 2 ) is an attempt to include the tower of higher local condensates of (q(D 2 ) n q) type. The change from purely local 
approximation f(z 2 ) = 1 to nonlocal condensates (NLC) with a smooth function f(z 2 ) that rapidly decreases for large 
z 2 modifies the QCD sum rule results for the moments of the pion DA: they become smaller. In particular, the initial 
NLC calculation [lj] gave (£ 2 ) = 0.25, and the model DA proposed in Refs. [H,[l9[ is 

€ R (0 = - V^e , (13) 

TT 

which is wider than the asymptotic DA, but narrower than the flat DA. The NLC method was elaborated in later 
papers, see Ref. [20j for a review. The problem of the NLC approach is that while it attempts to model the towers 
of (q(D 2 ) n q) condensates, the towers of (q~G n q) condensates are neglected. Recently, Chernyak [2l| gave a specific 
example in which the two towers exactly cancel each other. This means that NLC results may underestimate the 
value of (£ 2 ), and cannot exclude a possibility of large (£ 2 ) > 1/3 values for the second moment and (£ 4 ) > 1/5 
values for the fourth moment of the pion DA. On the other hand, it is quite possible that NLC argumentation is not 
completely wrong, and CZ results overestimate the values of (£ 2 ) and (£ 4 ). In particular, recent lattice calculations 
0,11 give (£ 2 ) w 0.29 and 0.27, respectively, at the scale fi 2 = 4GeV 2 , which produces (£ 2 ) > 0.3 at scales ~ 1 GeV 2 , 
but definitely not (£ 2 ) > 0.4. 

A general comment is that converting the obtained values of (£ 2 ) and (£ 4 ) into models for the pion DA one should 
not restrict the models by the requirement that DA's must be given by a few first terms of the Gegenbauer expansion. 
There is no a priori principle justifying such a requirement: it is just an assumption which may or may not be true. 



III. PHOTON-PION TRANSITION FORM FACTOR 



The form factor i ? 7 * 7 * 7r o (q 2 , q 2 ) relating two (in general, virtual) photons with the lightest hadron, the pion, plays 
a special role in the studies of exclusive processes in quantum chromodynamics. When both photons are real, the 
form factor i ? 7 * 7 * T o(0, 0) determines the rate of the ir° — > 77 decay, and its value at this point is deeply related to 
the axial anomaly [22 ]. At large photon virtualities, this form factor has the simplest structure analogous to that of 
the form factors of deep inelastic scattering. As a result, comparing pQCD predictions [H, [24], [25], [26], [53, HI] with 
experimental data, one can get information about the shape of the pion distribution amplitude y> w (x). Experimentally, 
F 1 , 1 * 7I o(ci 2 , q 2 ) for small virtuality of one of the photons, q 2 ~ 0, was measured at e + e~ colliders by CELLO [291 ] . 
CLEO [3fJ and recently by BaBar [l2| collaborations. 



A. Perturbative QCD 

The behavior of photon-pion transition form factor at large photon virtualities was studied [HI, [24], [25| within 
perturbative QCD (pQCD) factorization approach for exclusive processes d i 11 [M|. Since only one hadron is 
involved, the 7*7*7r° form factor has the simplest structure for pQCD analysis, with the nonperturbative information 
about the pion accumulated in the pion distribution amplitude ip n (x). Another simplification is that the short-distance 
amplitude for 7*7* — > tt° transition is given, at the leading order, just by a single quark propagator. Theoretically, 
most clean situation is when both photon virtualities are large, but the experimental study of i^ 7 * 7 * w o(Q 2 , Q\) in this 
regime through the 7*7* — > ir a process is very difficult due to very small cross section. 

In the lowest order of perturbative QCD, the form factor for transition of two virtual photons with momenta q±, qi 
into a neutral pion with momentum p = q\ + qi is given by contribution of the handbag diagram (see Fig[2j) 

-,pQCD,LO/ 2 r ,2\ _ \/2 f 1 <Pir(%) 



Tf* ^ Q2 >- 3 J xq 2 + {l-x)q 2dX - (MJ 
Introducing the asymmetry parameter lo through q\ = — Q 2 (l + w)/2 and q\ = —Q 2 (l — ui)/2 gives 

^^-WJo l+^x-l) dx = W J(W) • (15) 
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FIG. 2: Handbag diagram for photon-pion transition form factor. 

Thus, if one would know the function J{oS), one could (in principle) obtain the pion DA ipn(x) by inverting the integral 
transform (|15[) . However, as already mentioned, this kinematics is very difficult for experimental study. If one of the 
photons is real, i.e. lo = 1, the leading-order pQCD prediction is 

Information about the shape of the pion wave function is now accumulated in the factor J. It equals 2 for the 
infinitely narrow ~ S(x - 1/2) DA, for asymptotic DA © we have J as = 3, while CZ model ([12} gives J cz = 5. The 
intermediate distribution (fT3|) produces J MR = 4. Thus, in addition to (£ 2 ), we have another measure of the width of 
the pion DA, the value of J. Note, that for the DA's listed above, the ordering in J values is the same as the ordering 
in (£ 2 ) values. However, the flat DA, for which (£ 2 ) is smaller than that for the CZ model DA, generates infinite 
value for J, which is a consequence of the fact that it does not vanish at x = 0. This divergence of the integral for 
J formally means that the standard perturbative QCD factorization approach is not applicable for the flat DA case. 
But, since the divergence is only logarithmic, one may hope that some minimal fix, like a cut-off, might be sufficient. 
The question, of course, is whether there is a real need to use the flat DA to describe the data on the photon-pion 
transition form factor. 



B. Logarithmic model 

Recent data on 7*7 — > it form factor reported by BaBar collaboration in Ref. [13] are well fitted by the formula 

Q 2 F r ,MQ 2 ) = V2f n (_!!_^°' 25 = ^Ml J^(Q 2 ) (17) 

for the range 4GeV 2 < Q 2 < 40GeV 2 . The most startling observation is that J cxp (Q 2 ) does not show a tendency to 
flatten to some particular value. The specific (Q 2 ) 13 power-law parametrization of the growth is, of course, a matter 
of choice. In this region, J cxp (Q 2 ) is in fact very close to the logarithmic function 

^(Q 2 )-ln(l + |^) , (18) 

if one takes M 2 = 0.6 GeV 2 , see Fig. [3] The two curves practically coincide for Q 2 > 15GeV 2 . 

It is easy to notice that J L (Q 2 ) can be obtained if one uses the flat DA (p n (x) = f n and changes xQ 2 — ► xQ 2 + M 2 
in the pQCD expression for the 7*7 — * ir° form factor: 

As discussed above, the idea to modify propagators 1/fc 2 — > l/(fc 2 + M 2 ) in integrals over the light-cone momentum 
fractions is rather old. The parameter M in such modifications is usually treated as the average transverse momentum 
of the propagating particle. However, the immediate observation is that the value M = 0.77 GeV is a little bit too 
large to be interpreted in such a way. Furthermore, the 1/xQ 2 —> 1/{xQ 2 + M 2 ) modification is equivalent to bringing 
in, before the integration over x, a tower of (M 2 /xQ 2 ) n power corrections, i.e., higher twists. But it is known [32| 
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FIG. 3: Comparison of the function J exp (Q 2 ) corresponding to the fit of BaBar data (blue online) and logarithmic model 
function J L (Q 2 ) (red online). The asymptotic pQCD prediction J as = 3 is also shown (green online). 



that the handbag diagram, because of its simple singularity structure, cannot generate an infinite tower of power 
corrections. Indeed, the propagator of a massless quark in the coordinate representation is ~ ^(z 2 )~ 2 . Expanding the 
matrix element of the bilocal operator 

(0^(0)75 M*)\P) = 6MU=o + z 2 U(zp)\ z 2 =0 + (z 2 f^(zp)\ z 2 =0 + ... , (20) 

we see that twist-6 and higher terms cancel the singularity of the propagator. Hence, there are just two terms in 
the OPE for the handbag contribution: twist-2 term that has 1/Q 2 behavior and twist-4 term (corresponding to the 
-075 ^D 2 4> operator on the light cone) that gives 1/Q 4 contribution. Operators with (D 2 ) n - 2 do not contribute, and 
so there is no infinite tower of (1/Q 2 )™ terms. 



C. Light-front formalism and Gaussian model 

To investigate a possible mechanism capable of generating a cut-off at small x, let us write the 7*771° form factor 
in the light-front formalism. The required expression was given in the classic paper [23j of Lepage and Brodsky on 
exclusive processes in QCD. Namely, the two-body (i.e., qq) contribution to the 7*771"° form factor is given by 

* nfSUfl") = ^J^j %ltty-I ^ k ^- < 21 > 

Here, q± is a two-dimensional vector in the transverse plane satisfying q\ = Q 2 , ej_ is a vector orthogonal to q± and 
also lying in the transverse plane [23|, and the cross denotes the vector product. It can be shown that for the wave 
functions of ^(x, k±) = ip(x, k\) type we have [HJ 

F q 1 o(Q 2 ) = 1=1 — 0/ 0(x, k , ) k±dk± . (22) 

Following j2f|, we take the Gaussian ansatz for the fc^-dependence of the light-front wave function, which we write 
in the form 



tO/ , i 47T 2 ^(a;) / k\ \ 
* G (^ ; fej.) = J exp 

xxaV 6 V 2crxx J 



(23) 



where a is the width parameter and y w (x) is the desired pion distribution amplitude. The result for the form factor 
is then given by 



y/2 f 1 (f n (x) 



1 — exp 



xQ 2 y 

2xa ) 



dx . (24) 



It contains the 1/xQ 2 pQCD contribution and a correction term which makes the integral convergent in the region of 
small x. An important observation is that the correction term in the integrand of Eq. (|24p reflects the kj_ dependence of 
the nonperturbative pion wave function. In the Gaussian ansatz, this integrand term has an exponentially decreasing 
rather than a power behavior for large Q 2 . This fact alone is sufficient to assert that it cannot be classified as a 
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higher-twist term. It comes from contributions invisible in the operator product expansion, which only sees the terms 
that have a powerlike behavior in 1/Q 2 before integration over x. Representing this expression for the form factor as 

f^ (Q 2 ) = ^j g (qV), (25) 

we find that, for the flat DA ip v (x) — f n , the function J G (Q 2 , a) has the following large-Q 2 asymptotic behavior: 

jG(02,,T)=ln (?) +7£ + O((T/Q2) ' (26) 

where 7^ is the Euler-Mascheroni constant. Comparing this result with the function J L (Q 2 ,M 2 ) ([TO)) obtained 
through the M 2 modification of the pQCD 1/xQ 2 propagator, we conclude that they have the same (up to 0(1/ Q 2 ) 
terms) asymptotic behavior if 

M 2 

a=—e^, (27) 

which gives cr = 0.53 GeV 2 for M 2 = 0.6 GeV 2 . In fact, plotting J L (Q 2 , M 2 = 0.6 GeV 2 ) and J G (Q 2 ,a = 0.53 GeV 2 ) 
together, we observe that these two functions practically coincide in the whole region Q 2 > 1 GeV 2 we are interested 
in (see Fig. QJ. Comparison of the model curve with BaBar experimental data is shown in Fig. [5] 



J 
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FIG. 4: Comparison of the logarithmic model J L (Q 2 ,M 2 = 0.6GeV 2 ) (red online) and Gaussian model J G (Q 2 ,a = 0.53 GeV 2 ) 
(blue online). 



To check if the magnitude of a is in a physically reasonable range, let us calculate the average transverse momentum 
for this Gaussian model. We have 



(ki(x)) = J d 2 k ± k 2 1 ^(x,k 1 ) d 2 k ± dx^(x,k ± )j = 2ax(l-x) , (28) 

and, hence, 

= / {k 2 ± (x))dx = J . (29) 



Thus, yj (kj_) = 0.42 GeV, which is rather close to the folklore value of 300 MeV. One should also take into account 
that the wave function under consideration describes the valence two-quark Fock component of the pion, which is 
presumably smaller than other components. 

Thus, the magnitude of the M 2 -parameter of the logarithmic model is close to 3.3 (kj_) rather than to the value (kj_) 
expected from a naive substitution xQ 2 —> xQ 2 + k\ in the quark propagator. As we explained, such a change has no 
theoretical grounds in the case of the handbag diagram. The justification of the ad hoc modification xQ 2 — > xQ 2 + M 2 
used in our logarithmic model, as we have seen, is more complicated. 



D. One-loop pQCD corrections 



As already discussed, distribution amplitudes in general depend on the factorization scale /i, i.e. in principle one 
should always write: (p(x,fi). This dependence is induced by radiative corrections. The standard procedure in pQCD 



8 



Q 2 F 7r7r o(Q 2 ) (MeV) 




FIG. 5: Comparison of model curve (solid, red online) with BaBar experimental data. The asymptotic pQCD prediction 
Q 2 -F 77 *„.o(Q 2 ) = \/2fn is also shown (dashed, blue online). 



calculations involving pion DA is to start with an auxiliary quark-antiquark state in which the quarks are on shell 
and share the total momentum P in fractions xP and (1 — x)P according to the "bare" distribution amplitude 
<Po(x,m q ). Calculating radiative corrections for a specific process, e.g. for photon-pion transition form factor, one 
obtains logarithms ln((3 2 /m 2 ) accompanied by factors which may be converted into convolution of the lowest-order 
short-distance amplitude Tq(x) with the evolution kernel V(x,y) and the bare distribution amplitude (po(jj,m q ). 
Combining the evolution factor with bare DA, one obtains the expression in which Tq(x) is multiplied by "evolved" 
distribution amplitude ip(x,aQ), with a being some number, which is usually chosen in such a way as to minimize 
the size of that part of the corrections which was not absorbed into the renormalized (i.e. evolved) DA. One may 
also start with massless on-shell quarks, and use dimensional regularization to regularize mass singularities that result 
from ln(Q 2 /m 2 ) terms for m q = 0. Then the bare DA depends on the dimensional regularization scale and one 
gets ln(Q 2 /^ 2 ) evolution logarithms calculating corrections to the amplitude of the short-distance subprocess. 
The one-loop correction for the 7*7 — > 7r° form factor was calculated in Refs. (26l. I27I l28l| . with the result 



dx 



tp w (x) 



~x^ 2 ~\ 2^ 



•In x 



xhix 
2(1 -x) 
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In a^l In 



As advertised, the term containing the logarithm ln(Q 2 / fi 2 ) has the form of convolution 



U^A (30) 



1 r 
xQ 2 F 2n\2 



In x I = 



1 



(31) 



of the lowest-order term Tq(£, Q 2 ) — l/£,Q 2 and the kernel 



ITT 



X \ x — § 



9(£>x)[l + 



1 



1 - x ^ ' \ i- 

governing the evolution of the pion distribution amplitude. The "+" -operation is defined by 



(32) 



\F& x)]+ = Ffo x) - 6(£ - x) / F((, x) d( 



(33) 



When the probing momentum Q is much larger than the initial normalization scale /x, one deals with large logarithm 
ln(Q 2 /y« 2 ). The latter can be eliminated by taking [i — Q, and the expression is produced in which the evolved DA 
ip v (x,Q) is integrated with the remaining part of the correction. It is not guaranteed, however, that the resulting 
correction will be small, and the idea is to take fi = aQ with a chosen in such a way as maximally reduce the size of 
the a s correction. 

In the context of the present paper, we are interested in what happens when the bare DA is flat: ipo(x,n) = f w . 
Since in this case all integrals in (|30[) simply diverge, let us take a regularized version of the flat distribution amplitude, 
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namely the function 



with r being a very small parameter, say r < 0.1. Then Eq. (j30[) gives 



Jr(Q,(J>) 



1 



3tt 



- 9 + O(r) - - - 3 + — - r + O(r^) In 



Q 2 



(34) 



(35) 



It is clear that if we take /j, = Q, we will be left with a huge correction ~ (2a s /37r)/r 2 , i.e. ~ 60 (a s /7r) for r = 0.1. 
Since the coefficient in front of ln(Q 2 //i 2 ) is dominated by 2/r term, while the /i-independent piece is dominated by 
its 2/r 2 part, we can compensate the latter by taking ln(Q 2 //i 2 ) = 1/r. This corresponds to the choice 



/' 



)2 e -l/r 



(36) 



Thus, if we take r = 0.1 to model the flat DA, the optimal choice for fi is something like /i 2 = 10~ 4 <2 2 . Even for the 
highest Q 2 reached in BaBar experiment, this gives fi 2 = 0.004 GeV 2 , a scale corresponding to distances much larger 
than the pion size. Evidently, we cannot evolve the pion DA down to such small momentum scales. The evolution 
must stop at some /i 2 , ~ Aq CD . Thus, the flat pion DA becomes a DA at "low normalization point" n = /Xo ~ Aqcd, 
below which there is no evolution. Moreover, as we have seen in the example above, the radiative corrections do not 
induce visible 0(Q 2 ) additions to the renormalization parameter. Thus, in this "pQCD version" of the scenario with 
the flat DA, we deal simply with ip^fx). It does not evolve in the photon-pion transition amplitude, so there is no 
need to specify at which scale it is defined. 

Furthermore, writing the square-bracketed term in Eq. |35|) as [A(r) — B(r) ln(Q 2 /^ 2 )], we can fine-tune the 
coefficient a by taking a — exp[— A{r) / B{r)\ so as to completely eliminate the one-loop correction. Still, the resulting 
/j, 2 — aQ 2 will be very small, and there will be no evolution change in the shape of flat DA. In other words, in the 
pQCD version of flat DA scenario, there is no need to consider radiative corrections for the photon-pion transition 
form factor: they all are absorbed by the pion wave function. 

The photon-pion transition form factor was investigated in Ref. [H| using large- N c radial Regge model for resonances 
coupled to q\ and q2 photons. The results obtained in this way may be interpreted as a model with flat pion DA 
at low normalization point. In particular, ~ logQ 2 behavior was obtained for Q 2 F 1 - 1 ^o{Q 2 ) in the large-Q 2 region. 
The authors used the leading logarithm prescription 4>k{x) — * <t>Tt{%,Q) to make comparisons with experiment, i.e. 
the pion DA in their approach is not exactly flat for large Q 2 . 

When the pQCD version of the scenario with flat pion DA is applied to pion electromagnetic form factor, the analysis 
of radiative corrections is very similar. The conclusion is that there is no need to consider the one-gluon-exchange 
diagram: the gluon line should be absorbed into the soft wave function. After that, only the soft contribution remains, 
and the form factor should be calculated nonperturbatively. 




FIG. 6: Diagram for charmonium decay into two pions: the gluon lines cannot be absorbed into soft pion wave function 



This does not mean that the flat DA scenario excludes the diagrams with gluon exchanges for all processes. Consider 
charmonium decays into two pions, \ c — > 7T7t. The two gluons present in the lowest-order diagram cannot be absorbed 
into the pions' wave functions, so this diagram remains. In pQCD, it produces the same integrals of 0(£)/(l — £ 2 ) type 
that diverge for flat DA. Thus, one should write a more detailed expression involving the fc^-dependent light-front 
wave functions for both pions, which is a challenging problem for future studies. The description of charmonium 
decays is a well-known success of CZ approach: if one uses pion DAs close to the asymptotic one, the theoretical 
results are well below the experimental data. In case of unmodified propagators, the flat scenario gives divergent 
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results for these amplitudes, while propagator modification brings them down to finite values. It is interesting to 
check if the resulting values are close to CZ ones. 

One may argue that our logarithmic or Gaussian models for the lowest-order term are more or less equivalent to 
a simple cut-off of the s-integral at x = M 2 /Q 2 value, which is essentially larger than x ~ exp[— 1/r] values that 
are responsible for the dominant 1/r, 1/r 2 , 1/r 3 terms in the analysis above. If one simply imposes the cut-off at 
x = M 2 /Q 2 in the pQCD expression (|30|) . one would get powers of ln(Q 2 / M 2 ) instead of of powers of 1/r, and since 
ln(Q 2 /M 2 ) < 4 in our case, the asymptotically nonleading terms (especially (-9/2) contribution, see Eg. ([50)1 ) are 
essential. But it is not clear if a simple x cut-off in the pQCD expression is a correct prescription in the one-loop 
case. In particular, one may notice that the leading-order formula (|22[) of the light-front formalism can be formally 
written in terms of pion DA <p^(x,fjL) taken at [i = xQ (cf. Eq. Q): 

f%AQ 2 ) = x jf M*,/* = X Q) ■ (37) 

So, if the most important values are x ~ M 2 /Q 2 , then one should take fi ~ M 2 /Q (which is ia ~ Q e~^® 2 / M *\ 
compare with (|36p ). i.e. again a very small value for large Q. However, to check if this reasoning extends to the 
one-loop case, one needs to calculate one-loop corrections in the light-front formalism keeping the /c^-dependence in 
the perturbative part and then convoluting the result with fc^-dependent nonperturbative wave function(s), which is 
a task going well beyond the scope of the present paper. 



IV. SUMMARY 



In this paper we discussed a scenario in which pion distribution amplitude is treated as a constant for all values of 
the light-cone momentum fraction x. We indicated that several approaches, in particular QCD sum rules and lattice 
gauge calculations give the values for the second moment (£ 2 ) of the pion distribution amplitude that are compatible 
with this proposal. We emphasized that the standard practice of building the model pion DAs as a sum of two or three 
lowest terms of the Gegenbauer expansion is just an assumption. Such an assumption, however, excludes flat DAs from 
the start. We calculated the photon-pion transition form factor using the light-front formula of Lepage and Brodsky 
and incorporating a qq wave function that gives flat pion DA and has a rapid (exponential, for dcfiniteness) fall-off 
with respect to light-front energy combination kj_/x(l — x). We demonstrated that the use of such a wave function 
is numerically equivalent to 1/xQ 2 — > l/(xQ 2 + M 2 ) modification of the quark propagator, with the parameter M 2 
being more than three times larger compared to the average square of the valence quark transverse momentum. The 
characteristic feature of our result is logarithmic ~ In (Q 2 /M 2 + 1) growth with Q 2 of the combination Q 2 F j * 1 t T °(Q 2 )- 
Such a growth is indicated by recent data of BaBar collaboration fl2| . In this respect, it looks very important to 
check these results at other facilities, such as BELLE. 
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